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I MPACT PHENOMENA AS FAC T ORS IN THE EVOLUTION OF THE EARTH 

R. A. F. GRIEVE and E. M. PARMENTIER 

Department of Geological Sciences, Brown University, 

Providence, RI 02912, U.S.A. 

ABSTRACT 

Large terrestrial impacts result in appreciable ef- 
fects on the local geology and in regional seismic 
activity. Impact in recent geologic time, however, 
is n)t a significant global geologic process compared 
to many endogenic processes. This may not have been 
true in early earth history when the impact flux was 
considerably higher. It is estimated that 30-200 latge 
impact basins could have been formed on the early earth. 
These large impacts may have resulted in extensive vol- 
canisT. and enhanced endogenic geologic activity over 
large areas. Initial modelling of the thermal and sub- 
sidence history of large terrestrial basins indicates 
that they created geologic and thermal anomalies which 
lasted for geologically significant times. The role of 
large-scale impact in the biological evolution of the 
earth has been highlighted by the discovery of sidero- 
phile anomalies at the Cretaceous-Tertiary boundary and 
associated with North ‘merican microtektites. Although 
in neither case has an associated crater been identified, 
the observations are consistent with the deposition of 
projectile-contaminated high-speed ejecta from major 
impact events. Consideration of impact processes reveals 
a number of mechanisms bv which large-scale impact may 
indu.e extinctions. Although irpact has been largely 
ignored as a terrestrial geologic process, it is apparent 
that it must be considered in both the early geologic and 
more recent biologic evolution of the earth. 
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INTRODUCTION 

Data from bodies which have preserved portions of 
their earliest crusts, such as the moon, Mercury, Mars, 
and some of the Jovian and Saturnian satellites, indi- 
cate that impact cratering was a dominant process in 
early surface and upper crustal evolution. Although at 
a much diminished level, impact continues to influence 
the local geologic evolution of many bodies and, in 
particular, the surface evolution of bodies lacking an 
atmosphere. Ground truth data provided by terrestrial 
impact structures have played an important role in at- 
tempts to model impact phenomena. Geoscientists, how- 
ever, have largely dismissed or ignored impact as a 
geologic process on the earth. This is understandable 
as endogenic activity has done much to erase the ter- 
restrial record of impact. Nevertheless, approximately 
one hundred terrestrial hypervelocity impact structures, 
ranging up in diameter to 'v 140 km, are currently recog- 
nized. This paper discusses some possible implications 
of impact for the evolution of the earth, with particu- 
lar emphasis on early crustal and biological evolution. 

TERRESTRIAL IMPACT CRATERS 

Many of the details of the impact process can be found 
in papers in Roddy e£ al_. (1977), Shultz and Merrill 
(1981), and references therein and need not be repeated. 
Rocks subjected to hypervelocity impact and thus to 
shock pressures above their Hugionot Elastic Limit, 

'i 5-10 GPa for crystalline rocks, are characterized by 
shock metamorphic effects. Details of shock metamorphic 
effects can be found in French and Short (1968), 
Markushev (1981), and others. They are the principal 
criterion for identifying terrestrial impact structures. 
Relatively small terrestrial craters have a bowl- 
shaped form, with a breccia lens of allochthonous mate- 
rial below the floor of the crater and a depth/diameter 
ratio of the crater in the autochthonous target rocks 
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of 1/ 3—1/ 4 . The breccia lens, which may contain a basal 
pool of impact melt, is the result of the inward collapse 
of the initial or transient cavity formed by excavation 
and displacement of the target rocks by the cratering 
flow-field. At diameters > 2 km in sedimentary targets 
and > 4 km in crystalline targets, the final crater is 
characterized by depth/diameter ratios of < 1/10 and an 
uplifted central core of shocked rocks forming a ctn-ral 
peak and/or rings. The relatively shallow form is the 
result of uplift of the floor and collapse of the rim, 
which probably took place during and subsequent to tran- 
sient cavity formation. The mechanisms generally con- 
sidered responsible for modification are elastic rebound 
of the compressed floor and/or large-scale gravitational 
collapse of the transient cavity and surrounding area. 


Summaries of the geology and geophysics of many terres- 
trial craters can be found in Dence et^ aK (1977) and 
Masaitis et^ a_l. (1980). In some instances, the unusual 
geologic environment produced by large impact st ructures 
has influenced the deposition and exposure of ore depos- 
its; e.g., Ni-sulphides at the Sudbury structure (Morri- 
son, 1982) and uranium at the Carswell structure (Johns, 
1970). Significant reservoirs of hydrocarbons aie also 
located within some structures; e.g., Boltysh (Yurk e_t 
al . , 1975), Viewfie^d (Sawatzky, 1977), and Red Wing 
Creek structures (Brenan £t al_ . , 1975). Impact has a 
profound effect on local geology, disrupting the physi- 
cal and chemical equilibrium of the target rocks and 
resulting, in some cases, in structures with lateral 
dimensions greater than the largest volcanic constructs. 
For example, the 100 km Manicouagan structure, Canada, 
resulted from the release of 'V 10 2 J of energy virtu- 
ally instantaneously at one point on the earth's surface. 
This corresponds to 10 2 -10 3 times the energy released 
annually by all earthquakes. Material originally sev- 
eral kilometers deep was uplifted to »‘he surface, 'v 10 3 
km 3 of the target rocks were melted, and the effects of 
the impact are visible over an area of 2 x 10 1 * km 2 
(Grieve and Head, 1983). One regional effect of large- 
scale impact will be s°ismic activity. The fraction of 
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impact energy converted to seismic energy is 'v* 10 u and 
the relationship: log SE = 1.44 (M) + 5.24, where SE is 

seismic energy in Joules and M is surface wave earthquake 
magnitude (B3th, 1966), indicates that a Manicouagan- 
sized event would result in a magnitude 11 earthquake. 
Analyses of ground motions in nuclear explosions (Cooper, 
1977) suggest that peak displacements of a few centime- 
ters might be expected ^ 1000 km from such an impact 
event . 

IMPACT AND EARLY CRUSTAL EVOLUTION 

The crater production rate in the eirth-moor. system was 
extremely high and rapidly declining 4. 6-3. 9 by ago. 
During this early period, impact was a fundamental geo- 
logic process on the moon. Large impact basins form the 
basic topographic, tectonic, and stratigraphic framework 
of the ojon and impact was responsible for the character- 
istics c f the second order gravity field and upper crus- 
tal seismic properties of the moon (Grieve, 1980b). The 
overwhelming effect of impact on lunar rocks dating from 
this period is evidenced by ^ 90% of the returned lunar 
highlands samples being impact breccias or melt rocks, 
with Monte Carlo simulations indicating that the average 
lunar highland sample was affected by a minimum of 1-2 
impact events (Horz et al. , 1976). 

No known ter.estrial rocks date from the period of in- 
tense bombardment and the earth's crustal evolution dur- 
ing the first 800 my of its history is conjectural. The 
oldest terrestrial rocks have ages of 3. 5-3. 8 by, al- 
though ages of 4.1 by have recently been reported on re- 
worked zircons from Australia, ihe lack of an earlier 
crust may indicate that thermal and mechanical instabil- 
ities resulting from intense mantle convection and/or 
bombardment inhibited the preservation of continental 
. 6ialic crust. If a terrestrial crust and lithosphere 
did exist during a portion of the period of intense bom- 
bardment, the effects of impact might be, to a first 
approximation, similar to that on the lunar highlands 
crust. The number of impact craters produced can be 
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estimated from the cratering record in the lunar high- 
lands, corrected for the larger gravitational cross- 
section of the earth, the effect of the higher escape 
velocity on impact velocity, and the effect of higher 
gravity on final crater dimensions. 

Corrections for gravitational cross-section and impact 
velocity are functions of approach velocity. For ap- 
proach velocities of 6-10 tm s 1 , which is the model 
r.m.s. approach velocity fo k the residual planetesimal 
swarm (Wetherill, 1977), it is estimated that the earth 
will collect 2-4 times as many bodies of the same mass 
as the moon and the resultant craters will be 'V 1.5 
times larger, depending on the energy-diameter scaling 
relationship used. As the size-frequency distribution 
of primary craters approximates N “ D 2 , where N is a 
emulative number and D is diameter, this translates to 
a correction factor of 4.5-5 for the number of craters 
of equivalent diameter. The craterinq rate for the lu- 
nar highlands is ^ 6 x 10 3 km 2 for D > 4 km (Basaltic 
Volcanism Stud£ Project, 1981) and the lower factor of 
4.5 with N c D 2 gives 200 impact basins with D > 1000 
km on the earth in the period 4. 6-3. 9 by. This estimate 
is considerably higher than previous minimum estimates 
of 'v 3 n -60 basins with D > 1000 km. which were based on 
approach velocities of 15-20 km s (Frey, 1980; Grieve, 
1930a). This may be unrealistic as it requires a sig- 
nificant number of bodies in the residual planetesimal 
swarm, compared to the one or more of mass 10 i3 g needed 
to account for basin formation on the moon (Wetherill, 
1977). 

Given the number of impact structures, it is possible 
to estimate parameters such as the cumulative exogenic 
energy added to the early earth, the production of im- 
pact melt, and the extent of crustal fracturing. Even 
for minimal estimates of the number of impacts, the val- 
ues for additional exogenic energy and impact melt pro- 
duction are of the same magnitude as present-day proces- 
ses of internal heat losses and island-arc volcanism, 
respectively (Grieve, 1980a). Among the effects of 
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large-scale impact on an early terrestrial crust will be 
the formation of topography on the order of a few kilo- 
meters, enhancement of the sub-impact thermal gradients 
in the lithosphere and asthenosphere due to uplift of 
originally deep-seated materials, and the potential for 
subsequent eruption of basalt due to adiabatic decom- 
pression. Given that the earth has a relatively thin 
lithosphere, which is likely to have been even thinner 
during its earliest history, it is possible that large 
basin-sized impacts could bring the asthenosphere to the 
surface, leading to volcanism over large areas. Consid- 
eration of these effects led Frey (1980) to suggest that 
large impact basins ->n the early earth would produce 
proto-ocean basins floored by basaltic crust. 

Following impact, however, a number of processes act to 
modify large impact basins. These are: thermal contrac- 

tion and subsidence from the loss of post-shock and up- 
lift heat, topographic degradation by erosion and relax- 
ation, and loading of the basin. Some of these effects, 
as they apply to impact basins on the moon, have been ex- 
amined, with the result that lunar impact basins appear 
to have had appreciable effects, on timescales of 10 b -10 8 
y, on the thermal, volcanic, and tectonic history in and 
around the basins (Bratt et al., 1982; Solomon et al . , 
1982). Preliminary models of the evolution of a large 
terrestrial impact basin have been undertaken. A 1000 km 
basin, the form of which was constrained by data from 
smaller terrestrial structures (Grieve, 1980a), was 
formed in a model lithosphere 80 km thick with a thermal 
gradient of 20°C km 1 , overlying an asthenosphere with a 
thermal gradient of 0.5°C km 1 . Sub-impact temperatures 
were calculated as a function of time from finite differ- 
ence solutions of the transient heat conduction equation 
in a cylindrical geometry. The post-impact thermal anom- 
aly, created by post-shock waste heat and uplift of deep 
seated materials, decreases with time as heat is trans- 
ferred to the surface (Fig. 1). The initial thermal 
anomaly and the thermal anomaly after 24 x 10 6 y are shown 
in Figure 1. Substantial conductive cooling occurs on a 
time scale of 'v 10 7 years. 
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FIGURE 1. Vertical section of thermal anomaly beneath 
model 1000 km impact basin. R is radial distance and 
Z depth. Contours are for temperatures in excess of 
pre-innact thermal gradient. (a) Initial post-impact 
anomaly; (b) After 24 x 10 6 y. See text for details. 

Isostatic adjustment of the basin will occur on a time 
scale of 2y/p m gD (Parmentier and Head, 1981), where p^ 
and y are the density and viscosity of mantle material 
underlying the basin and g is the acceleration of grav- 
ity. For the model basin, isostatic rebound will occur 
in times on the order of 10 4 years for y=10 22 p.'ise. 

On longer time scales the basin floor will subside due 
to cooling and thermal contraction of underlying mate- 
rial. The thermal subsidence as a function of time, at 
several distances from the basin center, is shown in 
Figure 2. Thermal subsidence will be concentrated near 
the center, where a maximum subsidence of 2 km occurs 
after 'V 10 B years. In addition, the basin may be dy- 
namically loaded by volcanics, sediments and water, 
which will cause further subsidence. For example, a 
sediment load with density p $ s 2.7 g cm 3 and p = 

3.3 g cm 3 will cause a 1000 km basin with an original 
depth of ^ 3 km to isostatically subside by an addition- 
al factor of 1.8; i.e., the original floor will subside 
to a depth of 'v 5.4 km. This is in addition to thermal 
subsidence, which will contribute significantly to the 
thickness of sediments, particularly near the center. 
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•FIGURE 2. Thermal subsidence at various distances from 
the center of 1000 km basin, based on post-impact cool- 
ing model, as a function of time. 

The initial modelling described above does not account 
for several effects, which may be important in determin- 
* ing the magnitude and rate of thermal subsidence. First, 
loading of sediments into the basin will reduce the rate 
of heat loss to the surface and the rate of thermal sub- 
sidence. It has been suggested that if subsidence due 
to loading took place on time-scales shorter than the 
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decay of the enhanced thermal gradients below the basin, 
then the intra-basin materials: impact melt and breccia, 

post-impact volcanics, and sediments, may have been ther- 
mally metamorphosed and even partially melted to give 
sialic magmas (Grieve, 1980a). Second, the thermal anom- 
aly created by the impact will likely result in the up- 
ward transport of heat by solid-state convection and, 
therefore, the faster ceding of material at depth. Con- 
vective flow will also cause horizontal spreading of the 
thermal anomaly resulting in thermal subsidence less con- 
centrated near the basin center. Thus, convection will 
influence both the rate and radial distribution of ther- 
mal subsidence. These effects and the relative rates of 
subsidence and heat loss will be the subject of addition- 
al modelling. 

IMPACT AND RECENT EVOLUTION 

Recently, some authors have suggested relationships be- 
tween large impacts and such phenomena as magnetic rever- 
sals and certain plate movements (Glass et_ a_l. , 1979; 
Clube and Napier, 1982). These suggestions, however, are 
highly speculative. The most likely global effects of 
very large impacts will be coupled through the atmosphere 
and hydrosphere to climatic and biological changes. This 
thesis has been highlighted by the discovery of high sid- 
erophile abundances at Cretaceous-Tertiary boundary sites 
throughout the world (Alvarez et_ al_. , 19o0; Ganapathy, 
1980; and others) and apparently contemporaneous mass ex- 
tinctions. This geochemical evidence is consistent with 
the global dispersion of projectile-contaminated, high- 
speed, early-time ejecta (Grieve, 1982; O'Keefe and 
Ahrens, 1982) and it is difficult to argue against a ma- 
jor impact as the source of the siderophiles. Th« hy- 
pothesis of impact-induced extinctions is the sub.ect of 
discussion and is considered in detail in Geological So- 
ciety of America Special Paper 190 (1982). There is no 
known impact structure commensurate in size (^ 100-200km) 
and age (^ 65 my) with that required by the siderophile 
anomaly at the Cretaceous-Tertiary boundary. Three 
structures with D > 25 km in the USSR: Kamensk, Kara, 
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and Urst-Karska nave ages of 60 ± 5 my (Masaitis and 
Mashchak, 1982'' but they may not represent a significant 
increase over tne expected number of craters, provided 
they represent the total population of structures of this 
size and age (Grieve, 1982). The lack of an appropriate 
impact structure, however, may not be a problem, as Nd 
and Sr data suggest the presence of basaltic oceanic ma- 
terial within the boundary layer and thus an oceanic im- 
pact (Shaw and Wasserburg, 1982). 

A number of impact-induced extinction mechanisms have 
been suggested: (i) Injection of stratospheric dust, 

leading to the destruction of food chains (Alvarez et^ al , 
1980). The atmospheric flows generated by a large im- 
pact are sufficient to distribute stratospheric particles 
world-wide in a few hours (Jones and Kodis, 1982) and the 
debris cloud would reduce the amount of light below that 
required for photosynthesis for several months (Toon 
al. , 1982). (ii) Short-te,;m heating of the atmosphere 
and uppermost ocean layer by the kinetic energy of ejec- 
ta, followed by cooling due to dust (O’Keefe and Ahrens, 
1982). (iii) Poisoning by noxious chemicals (Hsi : , 1980) 
and/or alteration of atmospheric chemistry by lar^e quan- 
tities of NO and the rainout of HN0 3 (Lewis et_ al . , 

1982). It is clear that there are a number of impact- 
related effects which could sufficiently stress the bio- 
sphere to lead to mass extinctions. It is not yet clear, 
however, which if any was responsible for the Cretaceous- 
Tertiary boundary. 

Searches are cuirently underway for other siderophile 
anomalies in the geologic record. An Ir anomaly has been 
discovered in Antarctic ice and associated with the 1908 
Tungska event (Ganapathy, pers. comm.) and with North 
American microtektites , age 34 my (Alvarez et^ al_. , 1982; 
Ganapathy, 1982). This latter impact event has been 
. linked to the terminal Eocene biological and climatic 
changes (Alvarez et d., 1982) but, as with the 
Cretaceous-Tertiary, there is no known associated impact 
structure. The Popigai structure, USSR, is approximately 
the correct age, 39 i 9 my, and is sufficiently large, 
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D ^ 100 km (Masaitis et al., 1975), that it might be ex- 
pected to have produced global effects. Geochemical evi- 
dence, however, indicates that it is not the source cra- 
ter for the North American tektites (Shaw and Wasserburg, 
1982). 


Tab It 1^. Major extinction 

•venta and 

large impact atructurea 

of almllar age 

Extinction Event 

Age 

Impact Event 

Diameter 

tel 



Kara 

50 km 

60 t 5 my 

Sretaceoua-Tertlary 

65 1 1 my 

llrat-Karaka 

25 km 

60 1 5 my 



Kamenak 

25 km 

65 my 

Norian (L. Trlaaalc) 

214 t 6 my 

Manlcouagan 

100 km 

210 i 4 my 

Pernlan-Trlakilc 

245 1 5 my 






Charlevoix 

54 km 

360 4 25.0 my 

Fraanlan (L. Devonian) 

365 t 10 my 






Sll Jan 

52 km 

368 i 1.0 my 

Aahgllllan (L. Ordovician) 

421 t 4 ay 





Data Sourcca: Crleva (1982), Maaaltla at al.. (1980) and Sepkoakl (1982). 

Sepkoski (1982) recognizes five major, and a number of 
lesser, extinctions during the Phanerozoic. These are 
listed in Table 1, along with known large impact struc- 
tures with similar ages. It is tempting to speculate on 
correlations between specific impact structures and ex- 
tinctions; e.g., McGhee (1982). Such correlations, how- 
ever, are premature. Most impacts occur in the oceans. 

If all extinction events were related to impact, only 
•\> 25% should correlate with known land structures. The 
cratering record is incc nplete and age estimates of both 
extinction and impacts have attached imprecisions (Table 
1). Furthermore, there are a number of large structures 
which are not temporally associated with major biological 
perturbations (Grieve, 1982). Extinctions are a complex 
problem. They may be multi- jcsal and there are ques- 
tions on their severity, age , and time-spans. To estab- 
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lish firmly a correlation between an extinction event and 
a known crater would require more than a temporal overlap 
(Table 1). Some geochemical evidence in the form of sim- 
ilar siderophile anomalies in the stratigraphic record 
and the impact melts within the structure would provide 
much needed additional confirmation. 

SUMMARY REMARKS 

Although individual impacts may have significant effects 
upon local geology and seismlcally affect large regions, 
it is apparent that the global geologic effects of impact 
do not compare with modern terrestrial geologic processes. 
This may not have been the case, however, during early 
earth history when the impact flux was considerably high- 
er. As no recognizable evidence of the large multi-ring 
basins, which should have been produced at this time, is 
preserved, it is difficult to assess the extent to which 
large-scale impact influenced the early crustal evolu- 
tion of the earth. It seems likely that large impacts 
would have intensified and localized endogenic activity 
such as volcanism, sedimentation, and tectonism. Whether 
or not this activity was sufficient to modify impact ba- 
sins into stable areas of new sialic crust is not clear 
and requires additional detailed modelling. At this time, 
however, large-scale impact cratering must be regarded as 
having a potentially important effect upon early crustal 
evolution. 

The case of large Impacts affecting the biological and 
climatic evolution in more recent time is relative^ con- 
vincing. The evidence for a major impact at the 
Cretaceous-Tertiary boundary is restricted to geochemis- 
try. The discovery, however, of a similar anomaly asso- 
ciated with North American tektites, which arc generally 
recognized as having an impact origin, strengthens the 
case for the Cretaceous-Tertiary impact. Large impacts 
have several mechanisms through which they may affect 
both climate and the biosphere but interactions in the 
a tmosphere-hydrosp v er e-biosphere system are complex and 
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not \.ell known and, at present, no one extinction mecha- 
nism is clearly advocated. Civen the potential of large 
impacts to disrupt the biosphere, it is difficult to un- 
derstand why virtually all large impact events do not 
have associated extinctions* . This may be a recognition 
problem due to the lack of a complete fossil record or 
may mean that the biosphere must be already stressed be- 
fore impact can produce widespread extinctions. Alterna- 
tively, the Cretaceous-Terr 4 ary impact may represent a 
singular event. Although the exact relationship of major 
Impact events to mass extinctions may not be clear, it is 
apparent that large sporadic impacts must be giver con- 
sideration as a viable mechanism for mass extinctions and 
climatic changes on earth. Whatever the case, the dis- 
coveries at the Cretaceous-Tertiary will be the focus of 
considerable scientific study and, if nothing else, have 
forced a wide audience of scientists to consider impact 
as a process affecting the earth's evolution. 
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